Quantum efficiency of charge qubit measurements using a single electron transistor 
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The quantum efficiency, which characterizes the quality of information gain against information 
loss, is an important figure of merit for any realistic quantum detectors in the gradual process 
of collapsing the state being measured. In this work we consider the problem of solid-state 
charge qubit measurements with a single-electron-transistor (SET). We analyze two models: one 
corresponds to a strong response SET, and the other is a tunable one in response strength. We find 
that the response strength would essentially bound the quantum efficiency, making the detector 
non-quantum-limited. Quantum limited measurements, however, can be achieved in the limits of 
strong response and asymmetric tunneling. The present study is also associated with appropriate 
justifications for the measurement and backaction-dephasing rates, which were usually evaluated in 
controversial methods. 
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In quantum mechanics the Copenghagen's postulate 
assumes that a quantum measurement would instanta- 
neously collapse the state being measured onto one of 
the eigcnstates of the observable. This is the concept 
of perfect projective measurement. In practice, however, 
any realistic quantum detectors cannot realize such type 
of measurements. Actually, a realistic quantum measure- 
ment is a process of gradually collapsing the state being 
measured. In this context, the quantum efficiency is an 
important figure of merit for a quantum detector. To be 
more specific, let us consider the measurement of a two- 
state (qubit) system. Assuming the qubit is in an idle 
state which is simply a superposition of the logic basis 
states, but experiences no rotational operation between 
them. Further, we focus on a quantum measurement us- 
ing the so-called quantum non-demolition (QND) detec- 
tor, which only dephases the quantum coherence defined 
by the superposition, but does not flip the basis states. 
This situation coincides with the fact that the measure- 
ment operator is commutable with the qubit Hamilto- 
nian, which is one of the major critcrions of the QND 
measurement in general [l|, [2| • While for more general 
case the QND measurements of a qubit were discussed in 
Refs. [3, y , we restrict us in the present work to a sim- 
pler case as mentioned above with, however, a particular 
interest in the quantum efficiency. 

Now, consider the qubit measurements with a QND 
detector. During the (gradual) collapse process, one can 
get the measurement result only after some time until the 
signal-to-noise ratio reaches unity, owing to the stochas- 
tic nature of the elementary events leading to the collapse 
(such as tunneling or excitations in the detector). This 
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consideration actually defines a measurement time (r m ). 
On the other hand, the qubit state being measured would 
be inevitably dephased because of the detector's backac- 
tion, from which we can define a dephasing rate (I\j). 
Then, the quantum efficiency of a realistic detector is de- 
fined by tj = l/(2FdT m ), which can also be related with 
an interpretation of information gain versus information 
loss [a-0|- We say that a quantum detector performs 
quantum limited measurements when r\ reaches unity. 

In this work we consider a concrete realization of qubit 
measurements, say, a solid-state charge qubit measured 
with a single electron transistor (SET) [a-lllj. When the 
tunnel coupling between the charge states is quenched, 
the measurement falls into the scenario of QND type. 
Moreover, it is believed that the SET detector holds 
promisin g ap plications in solid-state quantum computa- 
tion [111 Il2| , whose measurement properties have been 
therefore received considerable attention in the past years 
[13l - ll9l |. While in the higher-order cotunneling regime 
the measurement can reach the quantum limit [r\ = 1) 
in principle [lot Il6j , it was found that the quantum ef- 
ficiency of the SET detector is rather poor in the weak 
response and sequential tunneling regime [12h14J . 

More recently, however, it was found that the signal- 
to-noise ratio in the power spectrum of qubit oscilla- 
tion measurements using the SET, another important 
figure of merit, can reach and even exceed the ideal 
value of quantum limited linear-response detectors [19(, 
i.e., the Korotkov-Averin bound [2Q|. To achieve such 
a result, the necessary conditions for the SET detector 
are an asymmetric tunnel coupling with the leads and a 
strong- response to the qubit. This result, together with 
some earlier investigations [171 13 ; provides a hint that 
the SET detector seems able to perform ideal quantum- 
limited measurement under similar conditions. In this 
paper, rather than the continuous measurement of qubit 



oscillations [13, LL9| 1 we w ^ u focus on the QND-type mea- 
surements for the idle (superposition) state of a qubit, 
to show how the quantum efficiency depends on the re- 
sponse strength and tunnel coupling to the leads and how 
a quantum-limited measurement can be achieved. The 
study will be associated with appropriate justifications 
for the measurement and backaction-dephasing rates, 
which were usually evaluated in controversial methods. 
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FIG. 1: Schematic setup of using single electron transistor to 
measure a solid-state qubit. Model (I) stems from the metallic 
SET, while model (II) can be realized with the semiconductor 
nanostructure of quantum dot SET. 



with different coupling strengths to the leads, i.e., ^l(r) 
and rw R y In model (II), the SET dot level is within 
the voltage window only for qubit in state \b). It lo- 
cates outside the voltage window if qubit in state \a), 
due to the stronger Coulomb interaction with the qubit. 
While model (II) is purely a strong response detector, 
we parameterize the response strength of model (I) as 
follows: Tl(T' l ) = (1 ± Ofi, T R (T' R ) = (1 ± Qf R , 
and 7 = Tr/Fl. Then, £ and £ characterize the re- 
sponse strengths of the detector to the qubit, while 
^L(R) = (^L(R) + r^m)/2 describe the average cou- 
plings. These rates, originally, are related with a num- 
ber of microscopic quantities, such as the tunnel-coupling 
amplitudes and the density-of-states of the leads through 
the Fermi's golden rule [T2|. In particular, the above 
identification of the rates is also associated with the con- 
sideration of large-bias voltage and low temperatures, 
which make the Fermi function be unity. 

Conditional Master Equation. — Following Ref. [19( 
and referring to the Appendix A of the present work for 
more details, a number-conditioned master equation can 
be formally expressed as 



Model Description. — Consider a charge qubit, say, an 
electron in a pair of coupled quantum dots with degen- 
erate dot levels E a = E^ and inter-dot coupling ampli- 
tude 0, measured by a single electron transistor. For the 



qubit, its Hamiltonian reads H, 
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fl(\a)(b\ + H.c). For the SET, in this work we consider 
two models [ljj, as schematically shown in Fig. 1. For 
both models, the qubit state can be discriminated from 
the different currents flowing through the SET detector, 
which are correlated with the qubit states |o) and \b) ow- 
ing to the Coulomb interaction H luii = Un c n a . Here, n c 
and n a are the number operators of the SET central dot 
and the qubit target state \a). 

In practice, there exist metallic and semiconduc- 
tor realizations for single electron transistors. For 
scmiconductor-quantum-dot-based SET, the models 
schematically shown in Fig. 1 can be understood straight- 
forwardly, by noting that the dot levels are discrete and 
assuming only one level in the voltage window. Accord- 
ingly, there are only two charge states, say, an empty and 
single-electron-occupation states. In the metallic case 
fl2| . the central dot is an island with dense energy levels. 
However, owing to the strong Coulomb blockade effect, 
it can still involve only two charge states in the transport 
process. Formally, these two charge states together with 
their couplings to the electrodes can be described simi- 
larly as the semiconductor quantum-dot SET. Particular 
treatments based on Fermi's golden rule for the metallic 
island-electrode coupling rates are referred to Ref. [12j . 

In our present work, we do not specify what type of 
SET is concerned, but instead present a unified study 
for two different models, i.e., (I) and (II) in Fig. 1. In 
model (I), the SET dot level is always within the volt- 
age window, regardless the qubit state in \b) or \a), but 
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Here, p^ describes the state of the system, say, the qubit 
plus the SET dot, conditioned on the electron number 
"n" tunnelled through the SET detector (i.e., the num- 
ber of electrons arrived to the right electrode). The Liou- 
villian £ is a commutator related to the system Hamilto- 
nian. The superoperators IZj describe electron transfer 
between the SET dot and leads, with explicit forms given 
in Ref. |19j . To implement this number-conditioned mas- 
ter equation in practice, we need to specify the state ba- 
sis. For the both models described above, we choose the 
same state basis: |1) = \0a), |2) = |06), |3) = |la), and 
1 4) = 1 16). In this notation |0(l)a(6)) means that the 
SET dot is empty (occupied) and the qubit is in state 
|a(6)). Explicitly, for model (I) we have [19| 
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while for model (II) it reads 
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In these equations, p< nR > means that the state is condi- 
tioned on the number of electrons tunneled through the 
right junction of the SET, which is more explicitly labeled 
here by "hr", instead of "n" in Eq. (p}. Also, differ- 
ing from Eq. ([!} which supports a bidirectional transport 
(j = +1 and j = -1), in Eqs. © and © only the for- 
ward unidirectional process (j = —1) survives. This is 
owing to the fact that the measurement is performed in 
a large-bias determined sequential tunneling regime. In 
this case, the SET-dot level that supports the transport 
current is deeply embedded between the Fermi levels of 
the two leads, leading thus to a strong suppression of 
the backward process. Finally, as explained in the part 
of introduction, we are interested in a QND-type mea- 
surement which allows us to set O = in Eqs. @ and 
©• 

Measurement and Dephasing Rates. — The measure- 
ment time r m , in realistic gradual collapse quantum mea- 
surement, is the average time needed to filter out the 
measurement signal from detector's noisy output. For 
model (I), it can be determined by a technique of wave- 
packet analysis |12| . which is actually equivalent to the 
counting-statistics technique employed in Ref. [6[ to ana- 
lyze the measurement rate with a strong response point- 
contact detector. Discrete-Fourier-transforming Eq. @ 
in terms of p(k,t) = Z™ =0 p M (t)e mk , with k e [0,2tt], 
yields 
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These equations can be split into three groups, i.e., Eqs. 
(4a)-(4b), Eqs. (4c)-(4d), and Eqs. (4e)-(4f). Consider- 
ing the characteristic solution proportional to e tult , for 
each group we can obtain two eigenvalues. Since for 
model (I) we are able to distinguish (read out) the qubit 
state only after relatively large number of electrons trans- 
mitted through the SET so that a measurement current 



is well defined, we need thus a solution only for small 
k, which dominantly contributes to the electron-number 
distribution function. Then, following Ref. [12], in the 
limit of k < 1, from Eqs. (4a)-(4b) and Eqs. (4c)-(4d), 
we can carry out the smallest two eigenvalues, formally 
expressed as uj^(k) = v^k + ^iv^f^k 2 , with p = a and b 
corresponding to the qubit states. Here, i> M are the wave- 
packet group velocities, being identical to the stationary 
currents J M associated with the qubit state \p), while / M 



The measurement time 
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are the respective Fano factors 

then reads [3, r m = (y/J K v° : + \/ f b v b 

For model (II), which represents a strong- response 
measurement, instead of the wave- packet analysis de- 
scribed above, a better way to determine r m is expected 
by the following consideration. Since for qubit in state 
I a) no electron can tunnel through the SET and arrives 
at the right lead, one then immediately knows the qubit 
in state \b) as soon as an electron is detected at the 
right lead. This is nothing but a single-shot measurement 
completed, and thus the measurement time is the aver- 
age waiting time for such a tunneling event. Applying 
the method of waiting-time analysis [22j , in Appendix A 
we present some details of calculating the measurement 
time, and obtain r m = 1/Tl + 1/Tfl. 

Now we discuss how to reliably determine the measure- 
ment backaction dephasing rate Td, for the present strong 
coupling models (I) and (II) (i.e. the large- U interaction 
between the SET and qubit). The underlying complexity 
originates from the fact that the off-diagonal element of 
pit), i.e., the simplest measure of qubit coherence, does 
not decay with single exponential rate in our case. Then, 
we may determine Yd differently as follows: 

(i) Regarding the slowest decay rate as the dephasing 
rate leads to Td = minpma;^], with uj^ the two eigen- 
values of Eqs. (4e)-(4f). Since this definition ignores 
the weight of each exponentially decaying component, it 
breaks down as the weight of the slowest decay compo- 
nent is the smallest one. 

(ii) Borrowing a technique in quantum optics, Yd can 
be defined better from the Glauber coherence function 
HI: g(r) = ([at(r)(7_(0) + h.c.]) M , where cr+ = \a)(b\, 
a^ir) is the operator in Heisenberg picture, and (• • -) ss 
means average with respect to the steady-state. Then, 
one defines T^ 1 = J °° g(r)dT. As we will see later, this 
definition suffers also to some extent the drawback of the 
method (i) above. 

(iii) In this work, in similar spirit of (ii), we propose 
to determine Y d simply as: T^ 1 = J °° p afc (T)dT/p ab (0). 
Differently, however, we assume and emphasize that at 
the beginning of measurement the SET dot is empty. We 
will see below that this identification is important, par- 
ticularly for small Yr, compared to the above Glauber 
coherence function method. 

Numerical Results and Discussions. — In Fig. 2 we 
plot T m , Yd, and the quantum efficiency r\ = l/(2IVr m ) 
against the setup asymmetry, respectively, for models (I) 
and (II). For the result of model (I), shown by the left- 
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FIG. 2: Measurement time r m , backaction dephasing rate 
Fd, and quantum efficiency r\ — l/(2TdTm), for models (I) 
and (II). For r m as shown in the upper panels, in addition to 
results from the wave-packet method for model (I) (solid line) 
and (II) (dotted line), the more desired result of model (II) 
is based on the waiting-time analysis (solid line). Plotted in 
the mid-panels are results of Td from, respectively, decay of 
pabij) with empty SET dot initially (solid line), decay of the 
Glauber correlation function (dashed line), and the smallest 
eigenvalue method (dotted line). Accordingly, the respective 
efficiency is plotted in the lower panels with the same type 
of lines. Parameters: Taking the SET dot level associated 
with the qubit state \b) as the reference of zero energy, and 
f_L and r_L as the energy units, we set [il = 50, jj,r = —20, 
and U — 40 for model (I) while U — 60 for model (II). And, 
for model (I), the response strength £ = £ = 0.9. 



column three panels of Fig. 2, the measurement time r m 
is carried out using the wave-packet analysis, while the 
dephasing rate T^ is computed respectively by the three 
methods (i)-(iii) outlined above. For Yd, we observe se- 
rious deviation from the desired result (solid line), indi- 
cating that either the smallest-eigenvalue method (dotted 
line) or the Glauber coherence function approach (dashed 
line) breaks down in this context, particularly for small 
Th/Pl- Since we are considering a strong coupling de- 
tector (i.e. with large U interaction) , the dephasing rate 
is dominantly determined by f £, being nearly indepen- 
dent of F ij/f l. Concerning the efficiency of model (I), 
we find that merely increasing the tunneling asymme- 
try Tji/Tl cannot reach the quantum limit in weak or 
even relatively strong response regimes. In addition to 
the tunneling asymmetry, the efficiency is also bounded 
by, and actually quite sensitive to, the response strength. 
This feature contradicts our intuition that an asymmetric 
SET is like a point-contact (PC) detector, by noting that 
an ideal PC detector has unit quantum efficiency which 
is independent of the response strength. We also notice 



that this conclusion differs from some statements in Ref. 
jig , but agrees well with the fact that the signal-to- noise 
ratio is considerably affected by the response strength 

MM- 

For the strong response detector model (II), we now 
demonstrate that the quantum limit can be achieved. 
The respective r m , T c i and r\ are plotted in the right col- 
umn of Fig 2. The solid- line in the upper panel plots the 
simple result r m = 1/Tl + \jY R from the waiting-time 
analysis, while the dotted-line shows the result from the 
wave-packet analysis as a comparison. For T^, shown 
in the mid-panel, we again find that the dephasing is 
dominantly caused by the tunneling events from the left 
lead of SET to its central dot, i.e., T d = T L /2. The 
dotted and dashed lines, respectively from the smallest 
eigenvalue and the Glauber coherence function methods, 
clearly show the extent of failure when applied to this 
model. Therefore, quite straightforwardly, the quantum 
efficiency is analytically obtained as: 



r\ = T R /(T L + T R ). 



(5) 



From this result, we see that the measurement is to be 
quantum limited with increasing T r /Yl 3> 1. 

Note that what this model (II) performs is actually a 
strong projective measurement since one can immediately 
conclude that the qubit is in one of the basis states if an 
electron tunnels through the detector. In terms of en- 
semble average times as used in the quantum efficiency 
definition (i.e., average over large number of transmit- 
ted and reflected events), the random tunneling through 
the left junction (equivalently, the random occupation of 
the SET dot) determines the dephasing rate of the qubit, 
while the detection in the right reservoir for the tunneled 
electron determines the measurement time. In this en- 
semble average sense, the qubit coherence is destroyed by 
the first jump, and the measurement rate is governed by 
the current r L r fl /(r L + T R ). 

For strong projective measurement, however, it may 
not be immediately clear how the fact whether the de- 
tector is quantum-limited or not would manifest itself in 
each single realization of measurement. We may explain 
this issue in terms of possible information loss in each sin- 
gle realization of projective measurement. Let us recon- 
sider the measurement process. That is, the left reservoir 
(source) emits a measuring electron, this electron first 
tunnels through the left junction, then through the right 
one, and is finally detected in the right reservoir (drain). 
From an information-theoretic point of view, the right- 
hand-side detector cannot distinguish whether there is a 
tunneling event from the left reservoir to the SET dot. 
This implies an information loss, and results in a non- 
unit quantum efficiency which makes sense even at the 
level of single projective measurement. To avoid this in- 
formation loss, one possible way is to detect the tunneling 
electron in the right reservoir with increase of the tunnel- 
ing asymmetry T r /Tl, as indicated by the above calcu- 
lation. Interestingly, based on the information-theoretic 
discussion, for model (II) with arbitrary T r /Tl, one can 



also reach the quantum limit by alternatively detecting 
electron tunneling through the left junction. This can be 
implemented by adding a nearby point-contact to detect 
the charge state of the SET dot, similar to the experiment 
by Gustavsson et al [2l|. In this way, the information- 
gain rate coincides with the dephasing rate, leading thus 
to a unit quantum efficiency. 

Conclusion. — To summarize, we have presented an 
analysis for the quantum efficiency of qubit measure- 
ments with a SET detector. Two models were analyzed: 
one corresponds to a strong response SET, and the other 
is a tunable one in response strength. We found that 
quantum limited measurements can be achieved only in 
a limiting case of strong response as schematically shown 
by model (II), together with an asymmetric tunneling 
setup. In the most range of response strengths, the SET 
detector cannot reach the quantum limit of efficiency. 
These results can be qualitatively understood by means 
of an information-theoretic interpretation, and are quan- 
titatively demonstrated with direct calculations by ap- 
propriately justifying the measurement and backaction 
dephasing rates. 
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Appendix A: Waiting-Time Calculation 

Conditioned on the number of electrons tunneled 
through the SET, rather than Eq. ([I]), we more explic- 
itly present the number-conditioned master equation as 



obtained and rewritten as 
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Here, the Liouvillian C is defined by £(•••) = [Hs, ■■■], 
where Hs is the system (i.e. qubit plus SET-dot) Hamil- 

tonian. We introduce operators 4 — C x (±£) a o 
with A = L, R and a\ (a c ) the creation (annihila- 
tion) operator of the SET central dot. The supcropcr- 

oters C x (=bC) are the generalized spectral functions: 
C[ ±} {±£) = j^ dtC x ±] (t)e ±iCt , where the bath corre- 
lation functions C x +) (t) = (fl x (t)f cX ) B and C x ~\t) = 

(fc\{t) fix) b ■ In these correlators, the average (• • -)b = 
Tr^[(' • • )pb] is over the local thermal equilibrium state 
(pb) of the SET leads, and /J A (f c \) are the leads creation 
(annihilation) operators from the tunneling Hamiltonian 
that describes the coupling to a c (a' c ). 

Notice that the usual unconditional state p(t) is re- 
lated to p^ n \t) simply as: p(t) = Xm/ 9 (*)• Based on 
Eq. (|A1[) . the master equation about p(i) can be readily 
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In large bias limit (i.e. sequential tunneling limit), the 
supcroperators are defined through 
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In these equations, A T( 
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is the Fermi function of the left (right) lead, while 
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Now we consider two counting schemes, and determine 
the respective waiting times until the first tunneling event 
happens. First, for a left-junction counting, conditioned 
on the result that no electron tunnels through it, the 
master equation reads 
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More explicitly, in the state basis |1) 
|3) = |lo), and |4) = 1 16) we have 
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For an initially empty SET dot, i.e., p(0) = |2)(2|, the 
solution reads pn(t) = p33(t) — pAi(t) — 0, and P22(t) = 
e Lt . Then, based on a quantum-jump concept we know 
the probability that the first electron tunnels through 

the left junction at time "£", p(t) oc Tr 
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r^e Tht . Accordingly, the waiting time for the first tun- 



neling event is 
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Second, for a right-junction counting and conditioned 
on the result that no electron tunnels through it, but 
regardless of the situation occurred at the left junction, 
the master equation reads 
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In the same state basis as above, we further have 
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Also, for initial condition p(0) 
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